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Abstract 

miRNAs are a class of non-coding endogenous small RNAs. They play vital roles in plant growth, development, and response 
to biotic and abiotic stress by negatively regulating genes. Mulberry trees are economically important species with multiple 
uses. However, to date, little is known about mulberry miRNAs and their target genes. In the present study, three small 
mulberry RNA libraries were constructed and sequenced using high-throughput sequencing technology. Results showed 85 
conserved miRNAs belonging to 31 miRNA families and 262 novel miRNAs at 371 loci. Quantitative real-time PCR (qRT-PCR) 
analysis confirmed the expression pattern of 9 conserved and 5 novel miRNAs in leaves, bark, and male flowers. A total of 
332 potential target genes were predicted to be associated with these 1 13 novel miRNAs. These results provide a basis for 
further understanding of mulberry miRNAs and the biological processes in which they are involved. 
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Introduction 

miRNA.s, which are found in animals and plants, are a class of 
19-24 nt non-coding small RNA molecules. They negatively 
regulate genes at the transcriptional and post-transcriptional level 
by cleaving the target mRNA and suppressing the translation of 
target mRNA [1]. miRNAs are encoded by MIR genes. In plants, 
MIR genes are transcribed into a pri-miRNA with a cap and a 
poly (A) tail and thereafter processed into a pre-miRNA, which is 
further cleaved into a mlRNA/miRNA* duplex. The last 
nucleotide of the 3' terminal in the duplex is methylated, and 
then the plant miRNA is loaded into the ARGONAUTEl 
(AGOl) complex leading to the cleavage of the target mRNA [1]. 
Increasing amounts of evidence have demonstrated that miRNAs 
play crucial roles in plant growth, development, and response to 
biotic and abiotic stress [2-4]. For example, miR156 regulates the 
transition of juvenile to adult m Arahidopsis thaliana [5]. miR172 
negatively regulates the cell fate specification in flower develop- 
ment of^. lhaliana [6]. miR160 controls the formation of root 
caps by targeting auxin response factors ARFIO and ARF16, both 
of which restrict the stem cell niche and promote columella cell 
differentiation [7]. 

Despite the importance of miRNA, the first miRNA, lin-4, was 
discovered untU 1993 from Caenorliabditis elegans [8]. In 2001, 
tens of miRNA were identified in several animal species by directly 
cloning and sequencing [9,10]. Since then, bioinformatic predic- 
tion and cloning have been used to identify many miRNAs in 
animals and plants [11-15]. However, although predictions that 
rely on the sequence can predict the conserved miRNAs easily, it 
difficult to identify species-specific miRNAs. The cloning method 
can only be used to identify small-scale miRNAs. In 2005, high- 
throughput sequencing technology was first used to sequence the 



small RNA libraries of A. thaliana and many miRNAs were 
identified [16]. This next-generation sequencing technology 
enables massive sequencing and detection of minimally abundant 
small RNA. It has become technique of choice for sequencing of 
the genome, transcriptome, and small RNA transcriptome [17- 
21]. A recent analysis of miRNA (based on Release 19, http:// 
www.mirbase.org/, August 2012) showed a total of 21,264 
miRNAs to be registered in miRBASE. This number is almost 
20 times than that in Release 6.0 (Release 6.0, ftp:/ /mirbase.org/ 
pub/mirbase/CURRENT/README, April 2005). The expan- 
sion of currently available miRNA information can be attributed 
to the development of this technique. 

Mulberry trees are widely planted in Europe, Africa, Asia, and 
the United States. This tree belongs to the genus Morus, family 
Moraceas, order Resales [22,23]. Mulberry leaves have been used 
to feed silkworms for silk production for about 5,000 years. In 
addition to being the sole nutritional source of the silkworm, 
mulberry tree have many other multiple uses. In particular, the 
secondary metabolites of mulberry plants are widely used as 
medicines [24—26]. In the current genomic era, the genome data 
are an important resource for gene identification and character- 
ization. The completion of the mulberry genome has allowed 
scholars to look at mulberry genes comprehensively [18]. Deep 
sequencing of transcripts can reveal many important gene 
products, such as miRNAs, which have been shown to be crucial 
to plant development and stress responses [2^]. For this reason, 
mulberry miRNAs were sequenced and analyzed in the present 
study. Three small RNA libraries of mulberry tissues (leaves, bark, 
and male flowers) were constructed and used for sequencing. 
Conserved, novel mulberry's miRNAs and their target genes were 
identified. The expression profiles of 9 conserved and 5 novel 
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mulberry miRNAs were confirmed in three tissues using stem-loop 
quantitative real-time PGR. These results expand our knowledge 
of the diversity and specificity of mulberry miRNAs and provide a 
basis for further understanding of the biological mechanisms that 
take place in mulberry plants. 

Materials and Methods 

Plant material and construction of small RNA library 

The wild mulberry species Moms notabilis grows in a pristine 
forest in Yaan, Sichuan Province, in southwest of China. This tree 
is located at 29°45.278' north latitude, 102°53.878' east longitude. 
It is a mull)crry species used for the genome sequencing. Young 
leaves, bark and male flowers from M. notabilis were collected in 
the spring of 2013. The samples were immediately frozen in liquid 
nitrogen and stored in — 80°C. No specific permissions were 
required for these activities. The field studies did not involve 
endangered or protected species. Small RNA libraries of M. 
notabilis were constructed as described elsewhere with a few 
modifications [27]. Briefly, the total RNA of three M. notabilis 
tissues was extracted using RNAiso plus (D9108A, Takara, China) 
in accordance with the manufacturers' instruction. Samples were 
then subjected to 15% denaturing polyacrylamide gel electropho- 
resis (PAGE). Small RNA fragments of 18-30 nt were separated, 
purified, and ligated with the 5', 3 'adapter sequentially. After 
reverse transcription and PCR, about 20 |ig products of three 
tissues were separately sequenced using lUumina HiSeq-2000 
(BGI-Shenzhen, China). 

Bioinformatics analysis of small RNA 

Raw data were filtered using a Perl script to delete low-quahty 
reads, chip adapter sequences, and contaminations. The sequences 
Sl8 nt of clean data were annotated in the Rfam database 
(Release 10.1) (http://www.sanger.ac.uk/software/Rfam) and 
Genbank non-coding RNA database (http://www.ncbi.nlm.nih. 
gov/) to remove non-coding RNA (rRNA, tRNA, snRNA, 
snoRNA) and degradation fragments of mRNA. The remaining 
sequences were aligned against miRNA database, miRBASE 
(Release 19) (http://www.mirbase.org/), and perfectiy matched 
sequences were considered conserved M. notabilis miRNAs. 

Prediction of novel miRNAs 

The unannotated small RNAs 18-30 nt in length were searched 
against M. notabilis genome. Novel M. notabilis miRNA were 
predicted using Mireap software (http:/ /sourceforge.net/ projects/ 
mireap/) with modifications basing on the default parameters: (1) 
the miRNA sequence length was 18-25 nt; (2) the miRNA 
reference sequence length was 20-23 nt; (3) the maximum copy 
number of miRNAs in any of the previous studies was 20; (4) the 
free energy of miRNA precursor was less than — 1 8 kcal/ mol; (5) 
the maximum space between miRNA and miRNA* was 300 nt; (6) 
the minimal base pairs of miRNA and miRNA* was 16; (7) the 
maximum bulge of miRNA and miRNA* was 4; (8) the maximum 
asymmetry of miRNA/ miRNA* duplex was 4; and (9) the flank 
sequence length of miRNA precursor was 20 nt. Hairpin 
structures of potential novel miRNA precursors were checked 
manually. The criteria were as follows: (1) the minimal folding free 
energy index (MFEI) of potential novel miRNA precursor was 
required be at least 0.85 [28]; (2) the asymmetric bulges of 
miRNA/miRNA* duplex were less than 3; (3) there were fewer 
than 4 mismatches between miRNA and miRNA*; (4) if the 
miRNA sequence did not fit these criteria, but corresponding 
miRNA* were detected, miRNA were also considered potentially 
novel miRNA. 



Expression of mulberry miRNAs and target genes as 
assessed using qRT-PCR 

Fourteen miRNAs were chosen for stem-loop RT-PCR in three 
mulberry tissues as previous described [29]. Briefly, each 1 (Xg total 
RNA was hybridized with a miRNA-specific stem-loop primer (10 
pmol). The hybridized miRNA molecules were reverse transcribed 
to cDNA in 10 |xL reaction using Reverse Transcriptase M-MLV 
(2641 A, TaKaRa, China) in accordance with the manufacturers' 
instruction. The resultant was then diluted three-fold and 1.5 |J,L 
cDNA was used as the template to perform the stem-loop RT- 
PCR with each miRNA specific forward primer and universal 
primer, as listed in Table S7. The reverse transcriptions for target 
genes were performed as follows. One microgram of total RNAs 
were reverse transcripted a 20 )xL reaction using PrimerScript RT 
Reagent kit with gDNA Eraser (RR047A, TaKaRa, China) basing 
on the handbook described. The resultant was then diluted four- 
fold and 1 nL cDNA was used as the template to perform the RT- 
PCR with each target gene primers, as Hsted in Table S8. The 
PCR reactions were performed in ABI Step One Plus (Applied 
Biosystems, USA) using SYBR Premix Ex Taq II (RR820A, 
TakaRa, China) as the following conditions: 95°C for 30 s, 40 
cycles of 95°C for 5 s and 60°C for 30 s. The 5.8S rRNA and 
ribosomal protein LI 5 gene were used as inner controls. AU 
reactions were assayed in triplicated. The relative expression level 
of miRNA was calculated using 2 method. 

Prediction of miRNA targets 

miRNA target prediction was performed by aligning miRNAs 
with the M. notabilis genes using a perl script, which was designed 
to predict the miRNA targets according to the criteria described 
by Allen and Schwab [30] [31]. These criteria were as follows: (1) 
There had to be no more than 4 mismatches in the miRNA/ target 
duplex (G— U pairs were considered 0.5 of a mismatch). (2) There 
had to be no more than 2 adjacent mismatches in the miRNA/ 
target duplex. (3) There had to be no adjacent mismatches at 
position 2-12 at the miRNA 5' terminal of the miRNA/ target 
duplex. (4) There had to be no mismatches in the position 10 and 
1 1 of the miRNA/ target duplex. (5) There had to be no more than 
2.5 mismatches in the position of 1-12 during the miRNA/ target 
duplex. (6) The minimal free energy of miRNA/ target duplex had 
to exceed 75% that of the miRNA when bound to its perfect 
complement. 

Results 

Small RNA in three mulberry tissues 

In order to identify the miRNAs in mulberr)' plants, total RNA 
(integrity ^8.0) was extracted from three mulberry tissues and 
three small RNA libraries were constructed for sequencing. These 
data have been deposited in NCBI/SRA database under accession 
number of SRP032829. 

A total of 11,752,747 reads (leaf), 11,491,921 reads (bark), and 
10,513,612 reads (male flower) were obtained by sequencing, as 
shown in Table 1 . After removing low quality sequences, adapters, 
and contaminated sequences, there were 10,992,174 (93.97%) 
clean reads £18 nt in size from leaves, 1 1,273,91 1 (98.55%) from 
bark, and 10,134,148 (96.83"'o) from male flowers. Among the 
clean reads £18 nt in size, the reads of miRNAs was 832,571 
(leaf), 1,130,016 (bark), 2,359,403 (male flower). The total reads of 
three tissues were subjected to analyze the size distribution as 
shown in Figure 1. The 21-24 nt small RNAs made up 77.87%, 
79.78%, and 81. 39%) of the reads from mulberry leaves, bark, and 
male flowers accounted, respectively. The most common size of 
small RNAs in leaves and bark was 24 nt, accounting for 38.75% 
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and 36.04% of the total, respectively. This was consistent with the 
size distribution patterns of small RNAs in Arachis hypogaea and 
Raphanus sativus [32,33]. However, the size distribution pattern 
of male mulberry flowers was different. The most common size for 
small RNAs in male flower tissue was 2 1 nt, which was consistent 
with Fragaria x ananassa and Pinus contorta [34,35]. 

Conserved miRNAs in mulberry plants 

In order to investigate the conserved miRNAs in mulberry 
plants, unique small RNAs from three mulberry tissues were 
aligned against miRNAs registered in the miRBase database 
(Release 19). Using the principle of sequence perfect matching, the 
present analysis identified 85 conserved miRNAs belonging to 3 1 
families, as shown in Table 2. There were 77, 70, and 70 miRNAs 
identified in leaf, bark, and male flower tissue, respectively. Among 
the 85 conserved miRNAs, 57 were common to all three tissue 
libraries (Figure SIB). The length distribution of 85 conserved 
miRNAs is shown in Figure SI A. A peak appeared at 21 nt 
(84.7%). 

In a broader evolutionary context, mulberry miRNAs were 
compared to those of seven other plants, including five dicotyle- 
dons {A. thaliana, Glycine max, Malus domestica, Populus 
trichocarpa, Ricinus communis) and two monocotyledons (Orym 
sativa and Zea mays). Of the 3 1 mulberry miRNA families, 24 
were conserved in the seven plant species. These miRNA were 
classified into well-conserved miRNA families. Prominent among 
them were mulberry miRlBOb, miR164a, miR167a, miR169a, 
miR390, and miR396b, which completely matched their counter- 
parties in the seven other plant species, suggesting those miRNAs 
were extremely conserved, and might play critical physiological 
roles in both dicotyledons and monocotyledons. However, 7 
miRNA families, miR482, miR529, miR858, miR4376, miR4414, 



miR4995, and miR5523, were found in only one or two plant 
species (Table 2). 

It has been reported that the sequencing frequency in lUumina 
technology was used to estimate the relative levels of expression of 
miRNAs [32]. The present data indicated that the conserved 
miRNA families were expressed across a vast range, from over 
10,000 reads to fewer than 10 reads in mulberry, as shown in 
Table 2. Of the 31 miRNA families, the reads of miR156, 
miR166, miR167, miR168, and niiR535 exceeded 10,000 in all 
three tissues. Eleven miRNA families (miR159, miRlGO, miR164, 
miR169, miR171, miR172, miR390, miR396, miR397, miR529 
and miR4376) had more than 1000 reads at least in one tissue. 
Seven miRNA families (miR162, miR393, miR395, miR398, 
miR399, miR408 and miR4414) had numbers of sequence reads 
ranging from 100-1000 at least in one tissue, and the remaining 
miRNA families (miR319, miR482, miR827, miR828, miR858, 
miR21 1 1, miR4995 and miR5523) had fewer than 100 reads in aU 
three tissues. 

Novel miRNAs in mulberry plants 

One of the greatest advantages of high-throughput sequencing 
is that this technology can be used to discover species-specific 
miRNAs. Here, mireap software was used with several criteria to 
identify the novel mulberry miRNAs, as described in method 
section. Using un-annotated sequences from mulberry three 
tissues, 262 novel miRNAs in 371 mulberry genome loci were 
identified, as shown in Table SI. Previous studies have reported 
that miRNA* sequence can be used as criteria for the identifica- 
tion of novel miRNAs. In the present study, 90 miRNA* 
sequences were discovered in the precursors of novel miRNAs. 
This suggested that those 90 novel miRNAs must be real mulberry 
miRNAs. Most of the 90 miRNAs* had much lower abundance 
than their partial strand miRNAs. However, several miR- 
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Figure 1 . Length distribution of small RNAs from mulberry leaves, bark, and male flowers. The X axis represents the length of small RNAs. 
The Y axis represents the relative frequency. 
doi:1 0.1 371/journal.pone.01 04409.g001 
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Table 1. Classification of smal 


RNAs In the 


mulberry leaves, 


bark, and male flowers. 








Category 


Leaves 




Bark 




Male flowers 






unique reads 


total reads 


unique reads 


total reads 


unique reads 


total reads 


Raw data 




11,752,747 




11,491,921 




10,513,612 


high quality 




11,697,570 




11,439,285 




1 0,465,855 


clean_reads [size > = 18 nt] 


3,383,828 


10,992,174 


3,848,791 


11,273,911 


2,851,747 


10,134,148 


match genome build2 


2,619,858 


9,394,886 


2,776,320 


9,122,054 


2,207,022 


8,601,051 


exon_antisense 


75,103 


1 95,404 


77,783 


178,514 


74,619 


1 60,433 


exon_sense 


175,011 


401,988 


151,119 


308,470 


150,431 


306,205 


intron_antisense 


140,785 


298,027 


153,833 


314,352 


108,403 


1 96,686 


intron_sense 


222,666 


865,955 


236,136 


714,362 


179,397 


511,639 


miRNAs 


26,956 


832,571 


31,463 


1,130,016 


25,805 


2,359,403 


rRNA 


114,561 


2,157,926 


162,589 


1,886,493 


111,055 


1,885,168 


repeat 


250,148 


525,805 


223,690 


414,606 


263,197 


435,246 


snRNA 


2,882 


13,281 


4,210 


17,357 


3,754 


18,088 


snoRNA 


1,058 


3,169 


1,148 


2,900 


831 


1,816 


tRNA 


12,054 


368,395 


20,001 


378,792 


13,571 


335,167 


unannotated small RNA 


2,362,604 


5,329,653 


2,786,819 


5,928,049 


1,920,684 


3,924,297 


doi:l 0.1 371/journal.pone.01 04409.1001 



NAs*(mno-miRnl66*, mno-miRn76a-2*, and mno-miRn82*) 
exhibited almost the same abundance with their partners in male 
flowers. Specifically, mno-miRn82* had higher abundance than its 
miRNA in leaf tissue. The 262 novel mature miRNAs and 371 
miRNA precursors were considered together, and the lengths of 
the miRNA precursors ranged from 65—357 nt with an average 
size of 142 nt. The majority were 65-189 nt in length, accounting 
for 87% of novel miRNA precursors. This was similar to those in 
A. hypogaea [32]. The minimal folding energy of the miRNA 

precursors varied from —225.7 21 kcal/mol with an average 

value of —56.50 kcal/mol, which was higher than that in tRNA (— 
27.5 kcal/mol) and rRNA (-33 kcal/mol) [36]. The majority of 
novel mulberry miRNAs identified in the present study were 
expressed at low levels with fewer than 100 reads. 

Identification of tissue-biased miRNAs 

Some miRNAs were expressed in a tissue-specific/biased 
manner. Documentation of such miRNAs has served as founda- 
tional information for functional studies. In the present study, 
tissue-biased miRNA families were investigated by normalizing the 
reads in different data sets. miRNAs were considered tissue-biased 
if there were twice as many normalized reads in one tissue than in 
the other two tissues. As shown in Figure 2 and Table S2, 5 leaf- 
biased miRNA families (miR4995, miR827, miR828, miR172, 
and miR390), 5 bark-biased miRNA families (miR535, miR396, 
miR393, miR395, and miR319), and 8 male flower-biased 
miRNA fiimilies (miR156, miR529, miR160, miR397, miR398, 
miR408, mLR169, and miR167) were identified in conserved 
mulberry miRNA families. Among mulberry novel miRNAs, 58 
leaf-biased, 84 bark-biased, and 72 male flower-biased miRNAs 
were identified, as shown in Table S3. Among these tissue-biased 
novel miRNAs, 8,6, and 1 5 novel miRNAs were only observed in 
leaf, bark, and male flower tissue, respectively (Table S4). These 29 
novel miRNAs had very few reads, ranging in number from 5- 
136. 



Detection of the expression of miRNAs using RT-PCR 

To confirm the expression pattern of the mulberry miRNAs, 9 
conserved and 5 novel miRNAs with different expression profiles 
were randomly selected for stem-loop RT-PCR analysis. As 
illustrated in Figure 3, miR827, miRn247, and miRnl84 were 
more abundantly expressed in leaves than in bark or male flowers. 
miR396 and miRn62 were highly expressed in bark. Five 
mulberry miRNAs, miR156, miR160, miR169, mLRn74, and 
miRnl88, were expressed predominantly in male flowers. The 
remaining 4 non-tissue-biased conserved miRNAs (miR159, 
miR162, miR164, and miR168) exhibited expression pattern 
nearly identical to the results of the analysis of sequencing data. 
Taken together, the results of stem-loop RT-PCR were consistent 
with the expression pattern of the tissue-biased miRNAs identified 
using high-throughput sequencing. 

Prediction of miRNA targets 

Plant miRNAs play important roles in biological processes by 
cleaving target mRNAs and suppressing the translation of target 
genes. In order to understand the biological functions of mulberry 
miRNAs, the target genes of 31 mulberry conserved miRNAs 
representing 31 families with high reads and 262 novel miRNAs 
were predicted using the methods described above. As fisted in 
Table S5 and Table S6, 89 target genes for 20 conserved miRNAs 
and 332 target genes for 1 1 3 novel miRNAs were annotated using 
the nr database. The majority of the target genes of conserved 
miRNAs were transcriptional factors, and as many miRNA targets 
were found to be conserved in mulberry plants as other plant 
species. These include miRNA-target pair associated with flower 
development, miR156-squamosa promoter-binding-like protein 
(SPL), miR159-MYB, miR166-homeobox-leucine zipper protein 
(HD-ZIP III), miR172-floral homeotic protein APETALA and 
pairs associated with root development, miR164-NAC domain- 
containing protein and miR167-auxin response factor 6 (ARF). 
Functional proteins were also identified as targets of conserved 
miRNA including mno-miR397 (laccase), mno-miR395 (sulfate 
transporter), and mno-miR390 (receptor-like protein kinase). The 
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Figure 2. Reads per million of conserved miRNAs in mulberry leaves, bark, and male flowers. A) miRNAs with >1000 reads; B) miRNAs 
with 30-1000 reads; C) miRNAs with <30 reads. The X axis represents different conserved miRNAs. The Y axis represents the reads per million for 
different miRNAs. RPIVI means reads per million. 
doi:1 0.1 371/journal.pone.01 04409.g002 



targets of novel miRNAs were mainly associated with protein- 
coding genes. For example, flavonol synthase/flavanone 3- 
hydroxylase, isoflavone 2'-hydroxylas, polyphenoloxidase, disease 
resistance protein, E3 ubiquitrn-protein ligase, basic proline-rich 
protein, aspartyl-tRNA synthetase, anthocyanidin 5,3-O-glucosyI- 
transferase, abscisic insensitive IB, chitinase-Iike protein and 
cysteine-rich receptor-like protein kinase. 

Detection of the expression of target genes using RT-PCR 

miRNAs would decrease the mRNA or protein level of their 
regulatrng target genes [1]. To explore whether the potential 
predicted miRNA targets could be regulated by miRNAs, the 
expression profiles of 10 target genes for one conserved miRNA 
and 5 novel miRNAs were investigated. As illustrated in Figure 3 
and Figure 4, the predicted target genes Morus015493 and 
Mo7'us018032 of miR156 were lowly expressed in male flower 
and highly expressed in leaf, which was opposite to the expression 
pattern of miR156. In addition, as shown in Figure 4 and Table 
SI, the predicted target genes Morus012124, Morus012122, 
Morus012121 of miRn51 also had the opposite expression pattern 
with miRn51. Four novel miRNA-target gene pairs 
(Morus008520 for miRn247, Morm019289 for miRn67, 
Morus011908 for miRn62, Morus002508 and Morus()14466 
for miRnl57) also possessed the opposite expression pattern with 
each other. 



Discussion 

MiRNAs are important components in regulating plant 
physiological processes [1]. In the past several years, abundant 
conserved miRNAs and species-specific miRNAs were identified 
by high-throughput sequencing because of its high-throughput 
capacity in the detection of large-scale miRNAs and high 
sensitivity in the detection of minimally expressed miRNAs. In 
this study, sequencing of the three mulberry small RNA libraries 
was performed using lUumina technology. After analyzing millions 
of small RNA reads from these RNA libraries, 85 conserved 
miRNAs belonging to 31 families and 262 novel miRNAs at 371 
loci were identified. After comparative analysis, several character- 
istics of conserved mulberry miRNAs were analyzed. A relation- 
ship between the degree of evolutionary conservation and the level 
of expression was observed in conserved miRNAs. As in previous 
studies, the majority of the conserved mulberry miRNA families 
identified here were evolutionarily conserved across plant species 
with high levels of expression. The less-conserved miRNA families 
(miR482, miR529, miR858, miR4376, miR4414, miR4995, and 
miR5523) showed lower abundance than the well-conserved 
miRNA families. The well- and less-conserved miRNA families 
may have evolved to play different roles in biological processes. 
The well-conserved mulberry miR164, miR167, miR156, 
miR172, miR159, miR166, miR171, miR172, and miR319 
targeted NACs, ARFs, SPLs, APETALAs, MYBs, HD-ZIPIII, 
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Figure 3. Relative expression of 14 mlRNAs in tKiree mulberry tissues (leaf, bark, and male flower) using stem-loop RT-PCR. The X axis 
represents different tissues. L indicates leaf tissue. B indicates bark tissue. MF indicates male flower tissue. The Y axis represents the relative 
expression level of miRNAs. 
doi:1 0.1 371 /journal.pone.01 04409.g003 



SCLs, AP2s, and TCPs, respectively. These transcription factors 
are very important to plant growth and development. For 
example, in^. thalinana, miR164 and miR167 afiect lateral root 
development and adventitious rooting, respectively [37,38]. 
miR167, miRi59, mlRlGO, and miR166 regulate the develop- 
ment of floral organs [39-46]. However, the targets of less- 
conserved miRNAs were mainly functional genes. It has recently 
been reported that miR482 and miR4376 target the NBS-LRR 
disease resistance gene and ACAIO, respectively, and that tliey 
played a role in disease resistance and reproductive growth 
[47,48]. Although well- and less-conserved miRNA families played 
different roles, both were found to be very important. Specifically, 
they cooperated to regulate tlie biological processes in plants. 

Many minimally expressed and species-specific miRNA have 
been discovered in the plant kingdom using liigh-throughput 
sequencing. This indicated that each plant has its own specific 
mlRNAs, which may play specific roles in physiological processes. 
In the present study, the characteristics of pre-miRNA including 
stem-loop structures and MFE, served as identification criteria for 
novel mulberry miRNAs. A total of 90 mulberiy miRNA* and 262 
novel mulberry miRNAs were identified. miRNA* is the product 



of dicer-like 1 (DCLl) and partial complementary to miRNA [49]. 
It was once considered to be degraded shortly after production and 
to have no roles in biological processes [50] . Recent studies have 
shown that miRNA* have important functions related to 
physiologically relevant levels, even though they are expressed at 
much lower levels than their miRNAs [51,52]. It is speculated that 
tlie 90 low abundance miRNA* identified here might play 
important roles in various biological processes in mulberry plants. 
The novel mulberry miRNAs exhibited features different from 
those of conserved mulberry miRNAs. The expression levels of 
most novel mulberry miRNAs were very low, and they mainly 
targeted functional proteins. The predicted potential targets of 
mulberry novel miRNAs were involved in cellular processes, 
metabolic processes, response to stimulus, and metabolism. 
Research into possible targets in other plants may provide 
important clues to facilitate understanding of the function of these 
novel miRNAs. In the present study, one target of mno-mLRn62 
may be the disease resistance protein. This suggests that this novel 
miRNA may play a role in disease resistance in mulberry plants. In 
mulberries, flavanone 3-hydroxylase (F3H) has been found to 
participate in the anthocyanin biosynthesis pathway. The possible 
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Figure 4. The expression profile of 10 target genes belonging to 6 miRNAs in three mulberry tissues (leaf, bark, and male flower) 
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expression level of target genes. MorusXXXXXX represented the target gene ID. 
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targets of four novel miRNAs (mno-miRnI4, mno-miRnl37, 
mno-miRnl76, and mno-miRn252) encoded F3H and may play a 
role in regulating anthocyanin biosynthesis. Further study into 
novel mulberr)' miRNAs may shed light upon their roles in 
mulberry biological processes. This may fill in the blanks with 
respect to current knowledge of biological processes involving 
conserved miRNAs. Novel miRNAs that work with the conserved 
miRNAs might regulate plant development and response to the 
environment more broadly and accurately than either set of 
miRNA alone. 

Conclusions 

This is the first comprehensive identification of conserved and 
novel miRNA in mulberry. The difierential expression of miRNAs 

and the prediction of their target genes provide a basis for further 
understanding of mulberry miRNAs and the biological processes 
in which they are involved. 

Supporting information 

Figure SI Length distribution (A) and tissue distribu- 
tion (B) of 85 conserved mulberry miRNAs. Numbers in B 
indicate the number of conserved miRNAs. 
(TIF) 

Table SI The information of 262 novel mulberry 
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Table S2 Tissue-biased conserved miRNAs and their 
target genes. 
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Table S3 Tissue-biased novel miRNAs and their target 
genes. 
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PCLSX) 

Table S5 Target genes of conserved mulberry miRNAs. 
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Table S6 Target genes of novel mulberry miRNAs. 
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